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We have performed a search for CP violation in a sample of -B^ — > DJ X decays corresponding 
to 5 fb~^ of proton-antiproton collisions collected by the DO detector in Run II at the Fermilab 
Tevatron Collider. New physics in B^ mixing could contribute a significant CP violating weak phase, 
which would be observed as a difference in the decay-time distribution for B^ — >■ B^ oscillated states 
versus that for — >■ B^. A fit to the decay-time distributions of the B^/B^ candidates yields the 
flavor-specific asymmetry a^^ = [— 1.7 ± 9.1(stat)^J^ g(syst)] x 10"'', which excludes CP violation 
due to new physics within the experimental sensitivity. 

PACS numbers: 11.30.Er, 13.20.He, 14.40.Nd 



I. INTRODUCTION 

The search for large CP violating (CPV) effects in the 
— i?j system is of special interest since their obser- 
vation would be a direct indication of new physics. A 



non-zero CPV weak phase (j)s arises from the phase differ- 
ence between the absorptive and dispersive parts of the 

- 5° mixing amplitude: 0^ = arg(-Mf2/rf2) IH, 
where and are the off-diagonal elements of the 

— i?° mass and decay matrices, respectively. A global 
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fit to various measurements interpreted in the context of 
the standard model (SM) yields the prediction (/)f = 
(4.2 ± 1.4) X 10-3 [ifl]. However, new physics, such as 
the existence of a fourth generation [llj , could contribute 
an additive phase such that (j>s = (/)f*^ + <j)f^- Re- 
cent measurements of cj)^^ in 5° — > J/ ipcl) decays by the 
CDF [l^ and DO [ll] collaborations differ from zero by 
approximately two standard deviations, which motivates 
further CP violation studies in decays. 

The CPV weak phase (f>s can be obtained from the 
flavor-specific asymmetry 

^B'>{t)^f +'^B°{t)^f 
c^ccoTdmg to a = tan (jjs , where AF^. and Am^ 
are the width and mass differences, respectively, between 
the heavy and light cigenstates of the mixed B'^ sys- 
tem. World average values of these quantities [l3| yield 
= {—SA'^l'^) X 10-3 dl]. Improved precision is 
needed to establish evidence of physics beyond the SM, 
which predicts a^^ = (0.0206 ± 0.0057) x 10-^ [l^. 

We present a measurement of a^^ using 5° — >• jX^Dj X 
decays (charge conjugate states are assumed throughout) 
reconstructed in proton- antiproton collisions collected by 
the DO detector between April 2002 and August 2008, 
corresponding to about 5 fb-^ of integrated luminos- 
ity [16]. This CPV study is complementary to those 
using inclusive dimuon events and B^ — >■ J/'ip'P decays. 
The time-integrated inclusive dimuon analysis [l7| does 
not distinguish between the various B hadrons and there- 
fore depends heavily on B'^ asymmetry results from the 
B factories and the determination of the &-quark frag- 
mentation fractions. In contrast, the present measure- 
ment allows a straightforward determination of the sam- 
ple composition, due to the partial reconstruction of the 
B° meson. The B° J/^cj} CPV measurements [H HI 
involve an analysis of the decay product transversity an- 
gles to separate the CP-even and CP-odd components. 
The present measurement does not require any angular 
analysis. Furthermore, it uses all the -B^ production and 
decay information available in an event, the former via 
initial-state flavor tagging, when possible, and the latter 
via an unbinned fit to the decay-time distribution. 

II. EVENT RECONSTRUCTION AND 
SELECTION 

The DO detector is described in detail elsewhere 
Charged particles are reconstructed using the central 
tracking system, which consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2 T superconducting solenoidal mag- 
net. An additional sing le- layer silicon microstrip de- 
tector called Layer [1^ installed immediately outside 
the beam pipe provides improved impact parameter res- 
olution and vertexing efficiency. Electrons are identi- 
fied by the preshower detector and liquid- argon/uranium 



calorimeter. Muons are identified by the muon system, 
which consists of a layer of tracking detectors and scin- 
tillation trigger counters in front of 1.8 T iron toroids, 
followed by two similar layers after the toroids '2^ . The 
solenoid and toroid polarities are reversed regularly, the 
latter allowing a determination of the muon charge asym- 
metries induced by the detector. The 5 fb-^ data sample 
used in this analysis is divided into two subsamples; the 
first 1.3 fb~^ is referred to as Run Ila and the remain- 
ing 3.7 fb^^ collected after the installation of the Layer 
detector is referred to as Run lib. 

Most of the sample was collected with single muon 
triggers. The reconstruction of the 5° — >■ n^DjX can- 
didates is as follows. The tracks were required to have 
signals in both the CFT and SMT. Muons were required 
to have measurements in at least two layers of the muon 
system. The muon segment was required to be matched 
to a track in the central tracking system and to have mo- 
mentum p(^+) > 3.0 GeV/c and pt(/^+) > 2.0 GeV/c, 
where pt is the momentum component transverse to the 
proton beam direction. 

All the tracks in each event were clustered into jets 
using the algorithm described in Ref. [2l| . The Dj can- 
didates were then formed from tracks found in the same 
jet as the muon candidate. Two fJ-^Dj^ final state sam- 
ples were reconstructed: /i+07r^ where 4> — > K^K^ and 
^i+K*°K- where K+tt-. The fi+cj^n^ recon- 

struction follows the technique described in Ref. [2^ . 
The (j) candidate was formed from two oppositely charged 
particles assigned the kaon mass (the DO detector is un- 
able to distinguish between kaons, pions and protons). 
The kaon candidates were required to have pt{K^) > 
0.7 GeV/c. The K~^K~ invariant mass distribution for 
the (p candidates in the final selected sample is shown in 
Fig.[T] The (j) candidate was required to have an invariant 
mass in the range 1.004 < M^K+K') < 1.034 GeV/c^ 
consistent with that of a </> meson. A pion candidate 
with charge opposite to that of the muon and pt{t^~) > 
0.5 GeV/c was then added to form the Dj^ meson candi- 
date. In the iX^K*^K~ decay mode, the D~ candidate 
was formed from three charged particles, one with the 
same charge as the muon and two with a charge oppo- 
site to that of the muon. The particle with the same 
charge as the muon was assigned the kaon mass and re- 
quired to have pt{K^) > 0.9 GeV/c. A more stringent 
requirement of pt{K~) > 1.8 GeV/c was imposed on 
the second kaon candidate to reduce combinatorial back- 
ground. The third particle was assigned the pion mass 
and required to have pt{tt^) > 0.5 GeV/c. The K'^7r~ 
invariant mass distribution for the K*^ candidates in 
the final selected sample is shown in Fig. [2l K*^ can- 
didates were required to have an invariant mass in the 
range 0.82 < M^X+tt") < 0.95 GeV/c^, consistent with 
the K*'^ mass. Details about the ii^K*^K~' analysis are 
available in Ref. [23| . 

The primary proton-antiproton interaction vertex was 
determined for each event. The average position of the 
collision point in the plane transverse to the beam was 
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FIG. 2: K'^TT invariant mass distribution for the fi^ K*°K 
sample. 



measured for each run and was included as a constraint. 
The precision of the primary vertex reconstruction for 
each event was on average about 20 /Ltm in the plane 
perpendicular to the beam direction and about 40 
along the beam direction. 

In both decay modes, a common Dj^ decay vertex was 
formed from the three Z?^ daughter tracks using the al- 
gorithm described in Ref. l24j]. To reduce combinato- 
rial background, the Dj vertex was required to have a 
displacement from the primary vertex in the transverse 
plane of at least four standard deviations. The cosine 
of the angle between the momentum and the direc- 
tion from the primary vertex to the Dj decay vertex 
was required to be greater than 0.9. The trajectories of 
the muon and Dj candidates were required to be consis- 
tent with originating from a common vertex (used as the 
decay vertex) and the fJ,'^D~ system was required to 
have an invariant mass between 2.6 and 5.4 GeV/c^, con- 
sistent with coming from a i?^ semileptonic decay. We 
define an angle between the combined fi^D~ momen- 



tum (an approximation of the momentum) and the 
direction from the primary vertex to the decay ver- 
tex. The cosine of this angle was required to be greater 
than 0.95 for candidates displaced from the primary 
vertex in the transverse plane by at least four standard 
deviations. These angular criteria ensure that the 
momentum is sufficiently aligned with that of its par- 
ent. The displacement and angular criteria give rise to a 
decay-time dependent reconstruction efficiency, which is 
discussed later. 

The B° selection was further improved using a like- 
lihood ratio method [1^ [2^ that combines a number 
of discriminating variables: the helicity angle between 
the and momenta in the (p or center-of- 
mass frame; the isolation of the fJ.'^Dj system, defined 
as / = p{fi~^ Dj ) / [p{fj,'^ Dj ) + 'Epi\, where the sum is over 
all tracks in the cone (A^)^ + (ATy)^ < 0.5 around the 
fj.^ Dj direction (0 is the azimuthal angle of the track, 
7/ = — ln[tan(6'/2)] is the pscudorapidity and 9 is the po- 
lar angle between the track momentum and the beam 
axis); the of the vertex; the invariant masses 
M{fj,+Dj), M{K+K-) {^x+(t>Ti- sample) or M{K+tt-) 
{li^K*'^K~ sample); and pt{K^ K') {ii^(f)-K~ sample) 
or pt{K~) {p^ K*^K~ sample). The final requirement 
on the likelihood ratio variable, ysei, was chosen to max- 
imize the predicted ratio S/^/S + B in a data subsample 
corresponding to 20% of the full data sample, where S 
is the number of signal events and B is the number of 
background events determined from signal and sideband 
regions of the A4{K'^K~tt~) distributions. Since j/sci is 
independent of the muon charge information, this op- 
timization does not influence the measured asymmetry 
and, therefore, the subsample was also used in the com- 
plete analysis. The numbers of Dj signal events, deter- 
mined from a fit to the K'^K^tt^ invariant mass distribu- 
tions (see Figs. El and m), are N{fi+(t)n-) = 81,394± 865 
and N{fi+K*°K~) = 33,557 ± 1,200, where the uncer- 
tainties are statistical. In approximately 1% of events, 
both a ^i^K*'^K^ candidate and a /i+(/)7r^ candidate 
were found. To avoid double counting, these events were 
removed from the ^'^K*^K~ sample. The events with 
two iX^(j)iT^ or ^'^K*^K~ candidates were removed from 
the analysis. 



III. FLAVOR TAGGING 

In order to measure the flavor-specific asymmetry a^^, 
it is necessary to distinguish between B^ — B^ and 
— > B^ oscillated states, which requires knowledge of 
the initial-state {production) and final-state {decay) fla- 
vors of the reconstructed B^ meson. The final-state b 
quark flavor is correlated with the charge of the muon 
in i?^ — fj,~^D~X semileptonic decays. The initial-state 



flavor, which provides additional information in the like- 
lihood fit that is used to extract a^^, was determined 
using an opposite-side tagging (OST) [13, algorithm. 
This algorithm relies in most cases on the reconstruction 



6 



o 

§"15000 
c 

0) 

> 

LU 

lOOOOt 



5000- 



L D0, 5 fb"' A 
~_ [I 71 sample / ■ 


-1- Data 
Pit 

\ - - - D3 ^ (\»K 

1 D' ^ 


" 1 1 1 1 ..r-'i ,...h.-V',L 1 


1 ■ ■ ■ ■ 



1.9 



2.1 



M,,JGevl^) 



FIG. 3: K TV invariant mass distribution for the jj.'^cjjTr 
sample with the solid line representing the mass fit result. 



25000 




FIG. 4: K+K'n 
fit result. 



M,^, (GeV/c=) 



invariant mass distribution for the 
sample with the solid line representing the mass 



of a second lepton (muon or electron) from the decay 
of the other b quark produced in the proton-antiproton 
interaction. This lepton appears on the side of the de- 
tector opposite to the reconstructed 5° meson (hence 
the term "opposite-side" tag) and its charge provides the 
flavor tag. When a second lepton cannot be identified, 
the OST algorithm attempts to reconstruct an opposite- 
side secondary vertex, in which case the charge is deter- 
mined from the tracks comprising the vertex. Only 21% 
of the events are tagged; the remaining 79% of events 
have neither a lepton nor a secondary vertex on the oppo- 
site side. Properties of the tagging lepton and secondary 
vertex tracks are incorporated in the tagging variable, 
dtag, which is assigned to each candidate. By defini- 
tion, the variable dtag is defined in the interval [-1,1]. An 
event with dtag > is tagged as an initial 5_quark and 
an event with dtag < is tagged as an initial b quark. A 
higher magnitude |rftag| corresponds to a higher tagging 
confidence. A sample of reconstructed fi^D*^ X 



decays was used to empirically determine the calibration 
function, I?((itag), called dilution. This function is used 
to calculate the probability Pcor = {p{dta.g) + l)/2 that a 
given 5° candidate has been tagged correctly. In events 
where no tagging information is available, the dilution is 
set to zero. 



IV. PROPER DECAY TIME 

The proper decay time of each i?^ candidate is de- 
rived from the measured displacement Lt of the B^ decay 
vertex from the primary vertex in the transverse plane. 
A Lorentz transformation of Lt into the B^ rest frame 
would yield the desired decay time. However, the unde- 
tected neutrino and other non-reconstructed particles in 
the semilcptonic B^ decay prevent the precise determi- 
nation of pT{B^g) needed to calculate the Lorentz boost 
factor. Instead, the combined transverse momentum of 
the n^Dj pair, pt{h^ D~), is used to calculate the visi- 
ble proper decay length (VPDL) 

I = M{B^) ■ [Lt ■ PT{^^^D-)]|[pT{^^+D-)]\ (2) 

where M{B^^) ^ 5.3663 GeV/c^ The proper de- 

cay length of each 5° meson is then ct{B'^) = IK,, 
where K, — pTilJ,'^ Dj) / pt{B1) is a correction factor that 
accounts for the missing momentum. Since K, is not 
known on an event-by-event basis, it was estimated from 
a Monte Carlo simulation, which included the PYTHIA 
generator [2^ interfaced with the evtgen decay pack- 
age [111 , followed by full GEANT [33| modeling of the de- 
tector response and event reconstruction. As large sam- 
ples were required to obtain sufficient statistical preci- 
sion, only generator-level information was used to de- 
termine the /C distributions. However, a sample of fully 
simulated events was used to verify that the difference be- 
tween generator-level K, distributions and those obtained 
using fully simulated/reconstructed events is negligible. 
A model of the muon trigger efficiency dependence on 
PxifJ-^) was included in the construction of the K. dis- 
tributions, which were obtained for each decay channel 
contributing to the signal sample. B^ semileptonic de- 
cays yielding an invariant ^jC^Dj mass that is close to the 
actual B^ mass have less missing momentum than those 
with lower M{fj,~^ Dj). Distributions of /C for a given de- 
cay channel in ten bins of M{iX^ D~) were used to exploit 
this fact, thereby reducing the uncertainty of the proper 
decay time associated with /C. 

The probability density function (PDF) for the B^ 
decay time is convoluted with the PDF describing the 
VPDL detector resolution and the PDF for the K, fac- 
tor. The decay-time PDF is then scaled by the re- 
construction efficiency, which was found for each decay 
channel using fully simulated events. In the ^"""^tt" sam- 
ple, the reconstruction efficiency for the decay channels 
contributing to the signal was then tuned to data by fix- 
ing the B1 lifetime, tbo — 1.470 ps [27|, and releasing 
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FIG. 5: The reconstruction efficiency as a function of VPDL 
for — > Dj X decays in tfie ^'^(f>Tr~ data samples. Plot 
(a) is for Run Ila and plot (b) is for Run lib. The width of the 
grey band indicates the range of the systematic uncertainty. 

the signal efficiency parameters in the decay-time fit to 
the data. In the ii'^K*°K~ sample, the overlap between 
Dj K*°K- and D" K*°n- candidates (see Fig.E]) 
prevented tuning to the data so the signal efficiency de- 
termined from the simulation was used as is, allowing 
tbo to float in the fit. For both decay modes, the re- 
construction efficiency for combinatorial background was 
determined from the data, as an adequate Monte Carlo 
model was not available. 

The decay-time dependent component of the signal effi- 
ciency was obtained by considering the events passing all 
the selection criteria except those dependent on the de- 
cay time. This efficiency is shown as a function of VPDL 
for the /i+f/jTT" and the iX^K*^K^ modes in Figs. [5] and 
[51 respectively. The VPDL range extends below zero to 
account for the VPDL resolution. The efficiency improve- 
ment visible in the low VPDL region in Run lib is due 
to the addition of the Layer detector. 

The functional form that best describes the efficiency 
in the positive VPDL region was chosen for the fit. An 
analytic function is required to properly normalize the 
decay time PDF. In the fj,~^K*'^K~ mode, where the 
signal efficiency was obtained from the simulation, the 




D0, Run Ila 

H K* K Monte Carlo 



-8.05 0.05 0.1 0.15 0.2 

VPDL (cm) 



- (b) 




D0, Run lib 




(I K* K Monte Carlo 







-8.05 0.05 0.1 0.15 0.2 



VPDL (cm) 

FIG. 6: The reconstruction efficiency as a function of VPDL 
for — >■ D~ X decays in the fi^K*'^K~ Monte Carlo sam- 
ples. Plot (a) is for Run Ila and plot (b) is for Run lib. The 
width of the grey band indicates the range of the systematic 
uncertainty. 

VPDL value for the minimum was allowed to float in the 
ffi. In the /i+(/)7r^ mode, where the signal efficiency was 
tuned to the data, it was necessary to fix the minimum 
efficiency point to be at VPDL= 0. A systematic un- 
certainty accounting for this was obtained by forcing a 
large negative variation in the efficiency in the negative 
VPDL region, as indicated by the width of the bands in 
Fig. [S] In the fj,^K*^K~ mode, the systematic uncer- 
tainty was obtained by simply varying the fit parameters 
within their uncertainties, as indicated by the width of 
the bands in Fig. [B] 

The VPDL uncertainty ai for a given event depends 
on the uncertainties of the production and decay vertex 
positions, which in turn depend on the track parameter 
uncertainties. The procedure described in Ref. [sij was 
followed to tune the tracking uncertainties using dijet 
and dimuon data samples. Realistic uncertainties were 
extracted from the track impact parameter pull distribu- 
tions, which were fitted by a single Gaussian function. 
The range for the fit was varied from between ztltr and 
±3cr. Parametrizations of the dependence of these un- 
certainties on the track momentum and polar angle were 
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thus obtained and applied to tune the track uncertain- 
ties in the analysis. The dimuon and dijet data samples 
yielded slightly different results. Using the tuned track 
parameter uncertainties, the vertex position uncertainty 
was also tuned, taking into account the tails, using the 
pull distributions of the vertex positions of J/?/' — > M^A*^ 
decays. Only the negative side of the pull distributions 
was fitted, as the positive side includes a long-lived com- 
ponent due to J/ 4! candidates from B decays. Separate 
vertex uncertainty tunings were obtained for data before 
and after the Layer installation. The tracking uncer- 
tainty differences observed in the dimuon and dijet sam- 
ples, combined with those associated with the Gaussian 
fit range variation, yield a 5% variation in the resulting 
VPDL uncertainty. A 5% systematic uncertainty on the 
VPDL was therefore included in the analysis. 



V. SAMPLE COMPOSITION 

There are two contributions to the Dj peaks in Figs. [3] 
and m the Bu,d,s — >■ f^^DjX decays (including muons 
originating from D(^s) and t+ decays) and the background 
occurring when the Dj^ meson originates from one b or 
c quark and the muon arises from the decay of another 
quark. The fraction of the peak that is background was 
determined from the data, as described later. The con- 
tributions from each B decay channel, accounting for 
branching fractions, were determined from the simula- 
tion. Measured values for the branching fractions were 
used. The exclusive branching fractions for scmileptonic 
decays to Dj, D^^ , and D*^ have not been 

measured. They were therefore calculated from the mea- 
sured branching fractions of the corresponding decays for 
-BJ] mesons, assuming the spectator model. 

The contributions to the /x+Dj signal prior to the ap- 
plication of the decay-time dependent criteria in each 
mode are shown in Table HI The relative contribution 
of each source varies with the invariant mass M{fi^D^) 
of the reconstructed n'^Dj system. For example, an 
event with a high M{fj,~^Dj) is more likely to originate 
from a direct — )■ ^'^v^D~ decay, rather than from a 
B^ ^^v^D*~ decay, where the intermediate D*~ de- 
cays to Dj X. The relative contributions for each source 
were therefore binned by invariant mass M[iX^ D'^) for 
an improved model of the sample composition. 

The mass PDF models the expected mass and width 
of the K^K^TT^ candidate for each source. Four sources 
were considered in the /x+^tt" sample (see Fig. [3]): the 
signal fj,^ Dj {—^ 4'''^~)] the accompanying mass peak 
due to ^^D^{-^ a- small reflection (less than 

1% and therefore not visible in the flgure) due to 
li^D~{^ A'+7r~7r~), where the kaon mass is misas- 
signed to one of the pious; and combinatorial back- 
ground. Five sources were considered in the iJ,~^K*'^K~ 
sample (see Fig. H]): the signal /i+£>7(-> K*°K~); the 
mass peak due to fi^D~{—i' K*'^K~); a reflection due to 



TABLE L Sample composition of the i-i'^Dg signal. 





Source 


H'^ (f>TT 


Fraction (%) 








20.69 




23.76 


B° 




63.26 




60.22 


B^s 




1.56 




1.65 


B° 




3.23 




3.16 


B° 


— > t'^D'Ut 


1.05 




0.25 




D-D+X 


0.68 




1.74 


B° 


DJDX 


0.68 




0.30 


B°s 


^ D+DX 


0.56 




0.30 


B+ 


D-DX 


3.57 




2.94 


B° 


D-DX 


4.72 




5.68 



a reflection due to /i+A^(^- K*^p)^ where the antiproton 
is mistaken for a kaon; and combinatorial background. 
The fractional contributions of these sources were deter- 
mined from the mass fits to the data. Note that the mass 
distributions shown in Figs. [3] and |4] are in fact averages 
of the mass PDFs for each event. 



VI. LIKELIHOOD FUNCTION 

AU events with 1.72 < M{(j)T:-) < 2.22 GeV/c^ 
(1.79 < M{K*°K-) < 2.25 GeV/c^) were used in an 
unbinned fitting procedure. While the signal events in 
both modes are confined to a narrower mass range, the 
events with masses outside the signal region are needed to 
accurately describe the VPDL distribution for the com- 
binatorial background under the signal peaks. The total 
likelihood C for N selected events is the product of like- 
lihoods Lj determined for each event j: 



N 



where 



[f^ 



(3) 



(4) 



The sum is taken over products of the probability den- 
sity functions for different sources of K~^K~'k~ candi- 
dates with fractions fi. The distribution of the VPDL 
uncertainty is described by Pp . The distribution of the 
likelihood ratio selection variable is given by P^'^K The 



mass PDF is given by P^ 



M{K+K^Tv~) 



Finally, Pf is 



fi+D-i^ K 



where the pion is mistaken for a kaon; 



the distribution of the tagging variable c?tag- These PDFs 
were determined from data. 

The PDF Pi for the measured VPDL / was constructed 
as follows. The formulae for the decay rates of neutral 
-B^ mesons were taken from Ref. [s^ assuming no direct 
CP violation (i.e., \Af \ — \Af\, where Af and Aj are 
the decay amplitudes). The asymmetry a^-^, modifies the 
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decay rates of mixed mesons as follows: 



TBOit)^f^Nf\Af\' 

[cosh {ATst/2) + cos(ATO5t)]/2, (5) 
rso(t)^/ = A^/l^/f (l-«/Je-i^=* 

[cosh {ATst/2) - cos{Amst)]/2, (6) 

[cosh {ATst/2) + cos{Amst)]/2, (7) 

[cosh {ATst/2) - coa{Amst)]/2, (8) 

where Nf is the normalization to the total number of 
mesons. Using the relationship ct{B'^) — IK, dis- 
cussed previously, these decay-rate PDFs can be written 
in terms of VPDL. All reconstructed events were divided 
into four samples corresponding to the final- and initial- 
state tags. The dilution of the initial-state tagging leads 
to a mixture of and i?° initial states in these samples, 
e.g., the sample tagged as 5° in the initial state and B'^ 
in the final state has the PDF 



'b<:'B° ^ ^ B'>{t)^fPcoi- 



+ r 



f (1 - Pc 



(9) 



where F 



s°(t)->/ 



f is the PDF in Eq. M and F^o 



the PDF in Eq. [T] As an example, the average VPDL 
PDFs for unmixed events (i.e., Pgogo and Pjjogo) and 
mixed events (i.e., P'^oga and PjjoBo) ^r the Run lib 
/i+07r^ sample are shown in Figs. [7] and [51 respectively. 
For the figures, the combinatorial background was sub- 
tracted using sidebands of the K~^K~t:~ invariant mass 
distributions. A contribution describing the background 
due to fake vertices around the primary vertex was in- 
cluded into the PDF. The slight differences in the PDFs 
for the mixed events in Fig. |8] give rise to a non-zero 
asymmetry a^^. 

The decay-rate PDF for the semileptonic decays B^ ^■ 
IJL'^D~X (the events comprising the peaks at 1.87 GeV/c^ 
in Figs. [3] and S]) is the same as for the B^ decays with 
the corresponding parameters changed accordingly. In 
particular, the B'^ semileptonic asymmetry a^^ was in- 
troduced and determined from the fit. 



VII. DETECTOR ASYMMETRIES 

The PDFs are modified to account for the detector 
charge asymmetries [J/7,, ^3] : 



=(1 + 7Adot)(l + <?^7^ro)(l + g7^fb)- 



(10) 



where /3 is the toroid polarity, 7 is the sign of the muon 
pseudorapidity, and q is the muon charge. Adet accounts 
for any asymmetry due to differences in the north {z < 0) 
and south {z > 0) ends of the detector (the proton beam 
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FIG. 7: The measured VPDL distribution for the Run lib 
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FIG. 8: The measured VPDL distribution for the Run lib 
/i^(^7r~ subsample tagged as B^ (B^) in the initial state and 
{B°) in the final state (i.e., mixed). The PDFs for B'^{t) -s- 
f and Bg{t) — >■ / are overlaid. 



travels from north to south). The range-out asymmetry 
Ay-o refiects the difference in acceptance of muons that 
bend towards as opposed to away from the beam line 
(for details, see Ref. [13 )• The forward-backward asym- 
metry Afb reflects the fact that positively charged muons 
tend to go in the direction of the proton beam whereas 
negatively charged muons typically go in the antiproton 
beam direction. Apj is a second-order correction to A^o 
that is non-zero only if Oy:^ and A^o are both non-zero. 
Aqp is a detector asymmetry between tracks bending to- 
wards 77 < and tracks bending towards 77 > 0. The 
likelihood for each event was multiplied by the detector 
charge asymmetry corrections for the muon from the B^ 
semileptonic decay and for the muon used for the OST 
when it is present. Both kaons in the (p-K^ signal sam- 
ple originate from the (f) decay and therefore have the 
same transverse momentum threshold. Since this is not 
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the case in the ^~^K*^K~ sample, the momentum depen- 
dent kaon reconstruction asymmetry [34] was taken into 
account. 



VIII. BACKGROUND DESCRIPTION 

The background where the Dj^ meson and the muon do 
not come from the same parent particle gives rise to fake 
vertices and a VPDL distribution that peaks around zero. 
Its shape was modeled by two Gaussian functions, and 
its contribution was estimated from decay-time fits to be 
approximately 8% for the total sample and 3% for the 
opposite-side tagged sub-sample for the ^i~^(/m~ mode. 
The corresponding contributions are 4% and 1% for the 
lj+K*'^K~ mode. The fake vertex background fraction 
is lower in the ^'^K*^K~ mode due to the higher kine- 
matic selection criteria. Several contributions to the com- 
binatorial background that have different VPDL distri- 
butions were considered. True prompt background was 
modeled with a Gaussian function. The long-lived com- 
binatorial background is dominated by misreconstructed 
heavy flavor decays. This background was modeled with 
an exponential function convoluted with the VPDL res- 
olution, including a component (w 60% of all long-lived 
background contributions) oscillating with a frequency 
of Arn^. An asymmetry parameter a^g was introduced 
into the PDFs for this component by analogy with Eqs.jS] 
and [71 This parameter absorbs possible asymmetries in 
the combinatorial background that were unaccounted for. 
The unbinned likelihood fit of the total sample was used 
to determine the various fractions of signal and back- 
grounds and the background VPDL parametrizations. 



IX. RESULTS 

The i?° oscillation frequency AttIs = 17.77 ± 
0.12 ps~^ [3^ was fixed in the fit. In the ii'^tpn" 
mode, the width difference was fixed to AP^ — 0.09 ± 
0.05 ps-i [13. In the n+K*^K- mode, the ratio 
ATs/Ts = 0.069t°;°^^ was fixed in the fit. First, the 
fi^(j)TT~ and fi^K*^K~ samples were fitted separately 
(see Table HI)) . The results for the charge asymmetries 
are consistent between these samples. The results of 
a fit to the combined /i+(/)7r~ and fj,~^K*'^K~ samples 
are shown in the last column of Table HT] The mea- 
sured value for the flavor-specific asymmetry is aj-^ — 
[— 1.7± 9.1] X 10^'^, where the uncertainty is statistical. 
The detector charge asymmetry A^o is observed to be 
non-zero, as expected [l3|- The other detector asymme- 
tries are consistent with zero, as is the physics asymme- 
try ttbg. There is no correlation between the detector 
charge asymmetries, either amongst themselves or with 
the physics asymmetries. There is a 5.2% (13.6%) corre- 
lation between Oy-^ and a^^ (ohg)- The non-zero value of 
aig is discussed in the next section. 



TABLE II: Asymmetries with statistical uncertainties. 











Combined 




X 10^ 


-7.0±9.9 


20.3±24.9 


-1.7±9.1 


d 

afs 


X 10^ 


-21.4±36.3 


50.1±19.5 


40.5±16.5 


abg 


X 10^ 


-2.2±10.6 


-0.1±13.5 


-3.1±8.3 


Afb 


X 10^ 


-1.8±1.5 


-2. Oil. 5 


-1.9±1.1 


Adct 


X 10^ 


3.2±1.5 


3.1±1.5 


3.1±1.1 


Aro 


X 10^ 


-36.7±1.5 


-30.2±1.5 


-33.3±1.1 


Apj 


X 10^ 


1.1±1.5 


0.2±1.5 


0.6±1.1 




X 10^ 


4.3±1.5 


2.0±1.5 


3.1±1.1 



X. SYSTEMATIC UNCERTAINTIES 

The sources contributing to the total a^^ systematic 

uncertainty of x 10~^ are hsted in Table [nil The 
largest contribution is due to the momentum-dependent 
kaon asymmetry. The fitted value of a'j^ is sensitive to 
this asymmetry due to the large D~ reflection under the 
Dj signal peak in the fi~^K*'^K~ mode. To account for 
this, and for the correlation between a'j^ and a^^, the 
kaon asymmetry was varied between zero and twice its 
nominal value and the likelihood fit was repeated in each 
case. Setting the kaon asymmetry to zero yielded a fitted 
value of — [—0.2 ± 16.5] x 10~^, whereas scaling it 
by two yielded aj^ = [81.2 ± 16.6] x 10"^^ The resuhing 
change in the fitted a^^ value is taken as the associated 
systematic uncertainty, given in Table IIIIl 

The same procedure was followed for the remaining 
sources of uncertainty. The signal fraction, that is, the 
number Dj and D~ signal events divided by the total 
number of events, obtained from the mass fits was varied 
up and down by Icr. The dilution of the opposite-side 
tagging algorithm was scaled by 0.9 and 1.1. The muon 
trigger efficiency was varied within its experimental un- 
certainties. The decay-time dependent reconstruction ef- 
ficiency as a function of VPDL was varied, as previously 
described. The VPDL resolution was scaled by 1.05 to ac- 
count for the tracking uncertainty differences previously 
discussed. The fake vertex background, where the Dj^ 
meson and the muon do not come from the same par- 
ent particle, was varied within its experimental range as 
was the prompt component of the combinatorial back- 
ground. The — i?g width difference APs and the os- 
cillation frequency Am^ were varied within their exper- 
imental uncertainties. The relative branching fractions 
(BP) for the exclusive semileptonic i?^ decays were var- 
ied within their predicted uncertainties in such a way as 
to keep the total inclusive fraction constant. To model 
a reduced signal fraction, the — > D^Df branching 

fraction was increased by la and the — >■ ^^Di*' X 
inclusive branching fraction was decreased by Icr. The 
total systematic uncertainty was determined by adding 
all the signed contributions in quadrature. 
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TABLE III: Systematic uncertainties. 



a{a%) X 10^ 


XVCHJll do V llllllC LI y O" Li HJ \J 


-1.24 


V\ Cii~\n Q^rTVi Tn P+ Q/^c*lp/n i^nr 
xVciUli db V i-l-li-l-lti Li y bCdicLl. UV ^ 


1 "in 


Signa,! frctction — Itr 


-u. / D 


SliP'Tial fraptinn -4-1 (T 

kJltilldl LL CLV^ L±Lj1± I _L Ly 


0.47 


Dilution scaled by 0.9 


-0.19 


Dilution scaled by 1.1 


0.21 


/i trigger efficiency low 


-0.03 


u trigger efficiency high 


0.00 


Decay-time dependent efficiency low 


0.15 


Decay-time dependent efficiency high 


-0.01 


VPDL resolution scaled by 1.05 


-0.03 


BF B° DjDt 


-0.00 




-0.10 


XVeiclLiVe JDI^ IDg 7 ^ U^LJg iUW 


n ni 

U.Ul 


Relative BF 5° — >• /i+i^,,_D7 high 


-0.05 


Fake vertex background low 


-0.13 


Fake vertex background high 


-0.04 


Prompt combinatorial background low 


0.01 


Prompt combinatorial background high 


-0.01 


AFs -Iff 


0.00 


AF, +la 


-0.01 


Arris — Icr 


-0.01 


Arris -l-lo" 


0.02 


Total 


+ i.4U 

-1.46 


XI. CONCLUSION 



In summary, using — > H^D^ X decays with 
D-^(t)n-, (j)^K+K- and D-^'k*^K-, 



j^*Q x^TT~, in combination with a decay-time 
analysis including initial-state flavor tagging, we mea- 
sure the asymmetry in mixed semileptonic decays 
to be a}^ = [-1.7 ± 9.1(stat)+J;5(syst)] x lO'^. This 
measurement supercedes the DO time-integrated anal- 
ysis of semileptonic B^ decays [sl], which yielded 
= [24.5±19.3(stat)±3.5(syst)] x lO'^. Our result is 
also consistent with the value a^-^ = (— 6.4± 10.1) x 10"'^ 
extracted [sB] from the DO time-integrated analysis 
of inclusive same-sign dimuon events [l7|. While the 
present result is the most precise measurement of the 
semileptonic asymmetry, improved precision is 

needed to establish evidence of CP violation due to new 
physics in B^ mixing. 
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